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Abstract 

We consider certain boundary conditions supporting soliton solutions in the generalized non-linear 
Schrodinger equation (AKNS r ) (r = 1,2). Using the dressing transformation (DT) method and the re- 
lated tau functions we study the AKNS r system for the vanishing, (constant) non-vanishing and the 
mixed boundary conditions, and their associated bright, dark, and bright-dark N-soliton solutions, re- 
spectively. Moreover, we introduce a modified DT related to the dressing group in order to consider the 
free field boundary condition and derive generalized N-dark-dark solitons. As a reduced submodel of 
the AKNS r system we study the properties of the focusing, defocusing and mixed focusing-defocusing 
versions of the so-called coupled non-linear Schrodinger equation (r—CNLS), which has recently been 
considered in many physical applications. We have shown that two— dark— dark— soliton bound states 
exist in the AKNS2 system, and three— and higher— dark— dark— soliton bound states can not exist. The 
AKNSr (r > 3) extension is briefly discussed in this approach. The properties and calculations of some 
matrix elements using level one vertex operators are outlined. 



Dedicated to the memory of S. S. Costa. 



1 Introduction 



Many soliton equations in 1 + 1 dimensions have integrable multi-component generalizations or more generally 
intcgrable matrix generalizations. It is well-known that certain coupled multi-field generalizations of the 
non-linear Schrodinger equation (CNLS) are integrable and possess soliton type solutions with rich physical 
properties (see e.g. [1, 2, 3, 4]). The model defined by two coupled NLS systems was earlier studied by 
Manakov [5]. Another remarkable example of a multi- field generalization of an integrable model is the 
so-called generalized sine-Gordon model which is integrable in some regions of its parameter space and 
it has many physical applications [6, 7]. The type of coupled NLS equations find applications in diverse 
areas of physics such as non-linear optics, optical communication, biophysics, multi-component Bose-Einsten 
condensate, etc (see e.g. [1, 2, 3, 8, 9]). The multi-soliton solutions of these systems have recently been 
considered using a variety of methods depending on the initial-boundary values imposed on the solution. 
For example, the direct methods, mostly the Hirota method, have been applied in [10, 3] and in ref. [11] 
a wide class of NLS models have been studied in the framework of the Darboux-dressing transformation. 
Recently, some properties have been investigated, such as the appearance of stationary bound states formed 
by two dark-dark solitons in the mixed-nonlinearity case (focusing and defocusing) of the 2— CNLS system 
[4]. The inverse scattering transform method (1ST) for the defocusing CNLS model with non- vanishing 
boundary conditions (NVBC) has been an open problem for over 30 years. The two-component case was 
solved in [12] and the multi-component model has very recently presented in [13]. The inverse scattering 
[12, 13] and Hirota's method [10] results on N-dark-dark solitons in the defocusing 2-CNLS model have 
presented only the degenerate case, i.e. the multi-solitons of the both components are proportional to each 
other and therefore are reducible to the dark solitons of the scalar NLS model. 

The r— CNLS model is related to the AKNS r system, which is a model with 2r dependent variables. As 
we will show below this system reduces to the r— CNLS model under a particular reality condition. Moreover, 
general multi-dark-dark soliton solutions of generalized AKNS type systems available in the literature, to our 
knowledge, are scarce. Recently, there appeared some reports on the general non-degenerate N-dark-dark 
solitons in the r— CNLS model [4, 14]. In [4] the model with defocusing and mixed nonlinearity is studied 
in the context of the KP-hierarchy reduction approach. The mixed focusing and defocusing nonlinearity 
2-CNLS model presents a two dark-dark soliton stationary bound state; whereas, three or more dark-dark 
solitons can not form bound states [4]. Here we show that similar phenomena are present in the AKNS2 
system. In [14] dark and bright multi-soliton solutions of the r— CNLS model are derived from multi-soliton 
solutions of the AKNS r system for arbitrary r in the framework of the algebro-geometric approach. 

In this paper we will consider the general constant (vanishing, nonvanishing and mixed vanishing- 
nonvanishing) boundary value problem for the AKNS r (r = 1,2) model and show that its particular com- 
plexified and reduced version incorporates the focusing and defocusing scalar NLS system for r — 1, and 
the focusing, defocusing and mixed nonlinearity versions of the 2-CNLS system, i.e. the Manakov model, 
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for r = 2, respectively. We will consider the dressing transformation (DT) method to solve integrable non- 
linear equations, which is based on the Lax pair formulation of the system. In this approach the so-called 
integrable highest weight representation of the underlying afhne Lie algebra is essential to find soliton type 
solutions (see [15] and references therein). According to the approach of [15] a common feature of integrable 
hierarchies presenting soliton solutions is the existence of some special "vacuum solutions" such that the 
Lax operators evaluated on them lie in some Abelian (up to the central term) subalgebra of the associated 
Kac-Moody algebra. The boundary conditions imposed to the system of equations must be related to the 
relevant vacuum connections lying in certain Abelian sub-algebra. The soliton type solutions are constructed 
out of those "vacuum solutions" through the so called soliton specialization in the context of the DT. These 
developments lead to a quite general definition of tau functions associated to the hierarchies, in terms of 
the so called "integrable highest weight representations" of the relevant Kac-Moody algebra. However, the 
free field boundary condition does not allow the vacuum connections to lie in an Abelian sub-algebra; so, 
we will introduce a modified DT in order to deal with this case. We consider a modified DT relying on the 
dressing group composition law of two successive DT's [29]. The corresponding tau functions will provide 
generalized dark-dark soliton solutions possessing additional parameters as compared to the ones obtained 
with constant boundary conditions. 

Here we adopt a hybrid of the DT and Hirota methods [7, 18] to obtain soliton solutions of the AKNS 
system. As pointed out in [17], the Hirota method presents some drawbacks in order to derive a general 
type of soliton solutions, noticeably it is not appropriate to handle vector NLS type models in a general 
form since the process relies on a insightful guess of the functional forms of each component, whereas in the 
group theoretical approach adopted here the dependence of each component on the generalized tau functions 
become dictated by the DT and the solutions in the orbit of certain vacuum solutions are constructed in a 
uniform way. The generalized tau-functions for the nonlinear systems are defined as an alternative set of 
variables corresponding to certain matrix elements evaluated in the integrable highest-weight representations 
of the underlying affine Kac-Moody algebra. In this way we overcome two difficulties of the methods. The 
Hirota method needs the tau-functions and an expansion for them, but it does not provide a recipe how to 
construct them. That is accurately solved through the DT method, which in its turn needs the evaluation 
of certain matrix elements in the vertex operator representation. We may avoid the cumbersome matrix 
element calculations through the Hirota expansion method. Since this method is recursive it allows a simple 
implementation on a computer program for algebraic manipulation like MATHEMATICA. 

The paper is organized as follows. In section 2 we present the generalized non-linear Schrodinger equation 
(AKNS2) associated to the homogeneous gradation of the Kac-Moody algebra sl^. The 2— component coupled 
non-linear Schrodinger equation (2— CNLS) is defined as a particular reduction by imposing certain conditions 
on the AKNS2 fields. In section 3 we review the theory of the DT, describe the various constant boundary 
conditions and define the tau-functions. In 3.1 we present the modified DT associated to the dressing group 
suitable for fee field boundary conditions. In section 4 we apply the DT method to the vanishing boundary 
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conditions and derive the bright solitons. In section 5 we consider the non-vanishing boundary conditions 
and derive the AKNS r (r = 1,2) dark solitons. Moreover, the general N-dark-dark solitons with free field 
NVBC of the AKNS2 arc derived through the modified DT method. In 5.3 we discuss the dark-dark soliton 
bound states. In 6 the dark-bright solitons arc derived associated to the mixed boundary conditions. In 
the section 7 we briefly discuss the AKNS r (r > 3) extension in the framework of the DT methods, and in 
the section 8 we discuss the main results of the paper. Finally, we have included the si (2) and sl(3) affinc 
Kac-Moody algebra properties, in appendices A and B, respectively. In C the matrix elements in sl(3) have 
been computed using vertex operator representations. 



2 The model 

The AKNS2 model can be constructed in the framework of the Zakharov-Shabat formalism using the hermi- 
tian symmetric spaces [19]. In [20] an affine Lie algebraic formulation based on the loop algebra g<E) C[A, A -1 ] 
of g € s/(3) has been presented. Here instead we construct the model associated to the full affinc Kac-Moody 
algebra Q = sl(3) with homogeneous gradation and a semi-simple element E"> (see appendix B). So, the 
connections are given by 

a = ^+E^4° ) +E*r^l+^, (2.1) 
i=i »=i 

i—1 i—1 i—1 i—1 

*+*-4°^ 2 - *+*r4tft + <t»c, (2.2) 

where , , and 2 are the fields of the model. In this construction we consider the fields as being 
real, however some reductions and complexification procedures will be performed below. Notice that the 
auxiliary fields and fa he in the direction of the affine Lie algebra central term C. The existence of this 
term plays an important role in the theory of the so-called integrable highest weight representations of the 
Kac-Moody algebra and they will be important below in order to find the soliton type solutions of the model. 
The connections (2.1)-(2.2) conveniently incorporate the terms with the constant parameters fli^ in order 
to take into account the various boundary conditions, so, these differ slightly from the ones in [20, 21]. The 
sl(3) Kac-Moody algebra conventions and notations are presented in appendix B. In the basis considered in 
this paper one has 

E® = \i2aH<f> (2.3) 

6 a=l 

The f3i in (2.1)-(2.2) are positive roots defined in (B.24) such that 

p\ = a\ +0:2; P2 = a 2, oi\,a2 = simple roots (2-4) 
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Notice that the connections lie in the subspaces defined in (B.12) 

A&g + gi, Beg + g 1 +g 2 , [D,g n ]=ng n (2.5) 

The zero curvature condition [<9 t —B ,d x — A] =0 supplied with the sZ(3) commutation relations 
(B.4)-(B.ll) and (B.25)-(B.31) provides the following system of equations 

2 



+ 



+^*+-2[2*t*T- l(ft.O)]*+ i = l,2 (2.6) 

2 

ft*r = -^*r + 2[2*+*7--(^)]*r, (2.7) 
d t 0i-d x 4> 2 = 0; (2.8) 

2 

= ^fiift, (ft.fi) = 20i + Q 2 , (y82.fi) = 2fi 2 + fii, fii, 2 = const .(2.9) 

i=l 

Notice that the auxiliary fields (/>i j2 completely decouple from the AKNS 2 fields ^f. The constant 
parameters Sli : 2 will be related below to certain boundary conditions and some trivial solutions of the 
system (2.6)-(2.7). 

The integrability of the system of equations (2.6)-(2.7), for fi 1)2 = 0, and its multi-Hamiltonian structure 
have been established [21]. A version of (2.6)-(2.7) for arbitrary r has recently been considered in [14]. The 
system of equations (2.6)-(2.7) supplied with a convenient complexification can be related to some versions 
of the so-called coupled non-linear Schrodinger equation (CNLS) [5, 1, 2, 3]. For example by making 

t -►-*«, [®f]* = Si $r = —fj, Si ipi, (2.10) 

where ★ means complex conjugation, fi e JR + ,Si — ±1, we may reduce the system (2.6)-(2.7) to the well 
known 2— coupled non-linear Schrodinger system (2 — CNLS) 

+ ^^l^f-^K^.fi)^ Vfe = 0, fc = l,2. (2.11) 

The term (/3fc.Q) is provided in (2.9) and the parameter fj, > represents the strength of nonlinearity 
and the coefficients 5j define the sign of the nonlinearity. The system (2.11) can be classified into three 
classes depending on the signs of the nonlinearity coefficients Si. For Si = 62 = 1 this system is the focusing 
Manakov model which supports bright-bright solitons [5]. For Si = S 2 = — 1, it is the defocusing Manakov 
model which supports bright-dark and dark-dark solitons [10, 12, 22]. In the cases Si = —S2 = ±1 one 
has the mixed focusing-defocusing nonlinearities. In this case, these equations support bright-bright solitons 
[23, 3], bright-dark solitons [24]. The defocusing and mixed nonlinearity cases have recently been considered 
in [4] through the KP-hierarchy reduction method and dark-dark solitons have been obtained. The system 
(2.11) has also been considered in the study of oscillations and interactions of dark and dark-bright solitons 
in Bose-Einstcin condensates [25, 26]. The multi-dark-dark solitons in the mixed non-linearity case are useful 
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for many physical applications such as nonlinear optics, water waves and Bose-Einstein condensates, where 
the generally coupled NLS equations often appear. 

The focusing CNLS system possesses a remarkable type of soliton solution undergoing a shape changing 
(inelastic) collision property due to intensity redistribution among its modes. In this context, it has been 
found a novel class of solutions called partially coherent solitons (PCS) which are of substantially variable 
shape, such that under collisions the profiles remain stationary [1, 27, 28]. Interestingly, the PCSs, namely, 
2-PCS, 3-PCS, r-PCS, are special cases of the well known 2-soliton, 3-soliton,..., r-soliton solutions of 
the 2-CNLS, 3-CNLS,..., r-CNLS equations, respectively [2, 3]. So, the understanding of the variable shape 
collisions and many other properties of these partially coherent solitons can be studied by providing the 
higher-order soliton solutions of the r-CNLS (r > 2) system considered as submodel of the relevant AKNS r . 
We believe that the group theoretical point of view of finding the analytical results for the general case of 
A-soliton interactions would facilitate the study of their properties; for example, the asymptotic behavior of 
trains of TV solitonlike pulses with approximately equal amplitudes and velocities, as studied in [16]. Notice 
that the set of solutions of the AKNS model (2.6)-(2.7) is much larger than the solutions of the CNLS 
system (2.11), since only the solutions of the former which satisfy the constraints (2.10) will be solutions 
of the CNLS model (2.11). This fact will be seen below in many instances; so, we believe that the AKNS r 
soliton properties with relevant boundary conditions deserve a further study. 

We will show that the three classes mentioned above, i.e. the focusing, defocusing and the mixed focusing- 
defocusing CNLS model can be related respectively to the vanishing, nonvanishing and mixed vanishing- 
nonvanishing boundary conditions of the AKNS r model in the framework of the DT approach adapted 
conveniently to each case. The (constant) non- vanishing boundary conditions require the extension of the DT 
method to incorporate non-zero constant vacuum solutions. Therefore, the vertex operators corresponding 
to the vanishing boundary case undergo a generalization in such a way that the nilpotency property, which 
is necessary in obtaining soliton solutions, should be maintained. Using the modified vertex operators we 
construct multi-soliton solutions in the cases of (constant) non- vanishing and mixed vanishing-nonvanishing 
boundary conditions. The free field NVBC requires a modification of the usual DT of [15] by considering 
two successive DT's in the context of the dressing group [29]. However, the same vertex operator generating 
the dark-dark solitons in the (constant) NVBC will be used in the free field NVBC case with a modified tau 
functions. 

3 Dressing transformations for AKNS2 

In this section we summarize the so-called DT procedure to find soliton solutions, which works by choosing 
a vacuum solution and then mapping it into a non-trivial solution, following the approach of [15]. For 
simplicity, we concentrate on a version of the DT suitable for vanishing, (constant) non- vanishing and mixed 
boundary conditions of the AKNS2 model (2.6)-(2.8). The free field boundary condition requires a modified 



5 



DT which is developed in subsection 3.1. So, let us consider 

, :h V"; -I'll' ^ pf,w 01,2 "►<>; pf L , pf R = const. (3.1) 

We may identify the NVBC (3.1) to certain classes of trivial vacuum solutions of the system (2.6)-(2.7): 

1) the trivial zero vacuum solution 

* J >ac = /^=0. i = 1 > 2 (3-2) 

2) the trivial constant vacuum solution 

*f, V ac = pf^0, .7 = 1,2. (3.3) 

Notice that (3.3) is a trivial constant vacuum solution of (2.6)-(2.7) provided that the expression Ej=i Pj Pj 
^(ft.O)] vanishes for any i = 1,2; which is achieved if Qi = fl 2 = ^, implying Y^j=i Pj Pj = 1^- 

3) the mixed constant-zero (zero-constant) vacuum solutions 

i)*t,vac = pt *0, *f <V ac = pt=^ (3-4) 

™)Kvac = pf =0, *£„ oc = P2 ?*0. (3.5) 

The first mixed trivial solution (3.4) requires 2pfpJ = 20i + O2, whereas for the second trivial solution 
(3.5) it must be 2p 2 p 2 = 2£1 2 + fii- 

The connections (2.1)-(2.2) for the above vacuum solution (3.3) take the form 

A- = E^ + piE^ + p- l E% + ptE^+p- 2 E% (3.6) 
73- ee £( 2 ) + pfE^ + p- lE % + p+E^ + p 2 E% ptp- 2 E^_ h ptp^Efl^ - 

(p+pj 1 )( J fff ) + iT<°>) (p+ P 2 V 2 )Hf\ (3.7) 

Notice that [A vac , B vac ] = and in order to define these vacuum connections it suffices to consider the 
constant values of PJ L / R \ in (3.1) related to one of the limits x — > ±00, say pj~ L = pj as above. These 
connections are related to the group element through 

A {vac) = a a .*(°)[*(°)]- 1 ; B(««=) =a t *(0)[*(0)]-l } ( 38 ) 

where 

*(0) =e xA««+tB»-«_ ( 39) 

The DT is implemented through two gauge transformations generated by 6± such that the nontrivial 
gauge connections in the vacuum orbit become [15] 

a = eiA^ie'i}- 1 + d x 0'l[6'l}- 1 (3.10) 
b = e'iB < - vac ^e'i}- 1 + dt&ii&ir 1 (3.11) 
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where 



9^ =exp(^a_J , 0+ = M~ 1 N ; M = cxp (<r ) , iV = exp ( ^ a„ j , 

\n>0 / \n>0 / 

where [-D,cr„] —na n . Therefore, from the relationships 

A = d x {V h )[* h }- 1 : B = dti^i^r 1 , * h = O^* (0) , 

one has 

where /i is some constant group element. 

One can relate the fields 'J* , (j>\ and </> 2 to some of the components in a n . One has 



(3.12) 



(3.13) 



(3.14) 



A 
B 



A vac + 



tr_i, E^ 1 ' + terms of negative grade. 



(3.15) 



+ 2 



C-l, 



<7-l,i5 (2) 



+ 



i=i 



+ 



(1) 



+ terms of negative grade 



(3.16) 



Taking into account the grading structure of the connection A in (2.5) we may write the a' n s in terms of 
the fields of the model. In order to match the zero grade terms of the both sides of the equation (3.15) one 
must have 

= -X>» + -ptH~ 1] +E(*r - Pi) E -S +E^i^- 1) - (3.i7) 

i— 1 i— 1 a— 1 

In the equation above, the explicit form of a a _ 1 in terms of the fields can be obtained by setting the sum 
of the (—1) grade terms to zero. Nevertheless, the form of the g°_ x will not be necessary for our purposes. 

Following the above procedure to match the gradations on both sides of eqs. (3.15)-(3.16) one notices 
that the cr_„ 's with n > 1 are used to cancel out the undesired components on the r.h.s. of the equations. 

From the equations (3.12)-(3.13) and (3.14) one has 



(A| M" 1 



A 



^(0)^(0)- 



A 



(3.18) 



where (A| and A \ are certain states annihilated by and Q >1 respectively. Defining 



°o = E < E ^ ] + E -o a ^l + E - a Hi°) + v c 

a>0 a>0 a=l 

and choosing a specific matrix element one gets a space time dependence for the field r\ 



-v 



(Ao| 



TO, 



^(0)/j^(0)-l 



|Ao) 



(3.19) 



(3.20) 
(3.21) 
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where we have defined the tau function tq and |Ao) is the highest weight defined in (B.15) . Next, we 
will write the fields tyf in terms of certain matrix elements. These functions will be represented as matrix 
elements in an appropriate representation of the affine Lie algebra sl(3). We proceed by writing the eq. 
(3.14) in the form 



exp - ]T <r_„ |A ) = \X ) r^ 1 



n>0 



(3.22) 



where the cqs. (3.12), (3.19) and (3.20)-(3.21) have been used. 

Then the terms with grade (-1) in the both sides of (3.22) can be written as 



-a-i \X a ) = 



t (x, t) 



(3.23) 



or equivalently 
/ 2 



(o)-il 



u |Ao> 



\i=l 



i=l 



a=l 



T (x, t) 



(3.24) 



Acting on the left in cq. (3.24) by E±p. and taking the relevant matrix element with the dual highest 
weight state < A Q | we may have 



Pi 



TO 



and 



i = 1,2. 



T 



where the tau functions , t are defined by 

T+( X ,t)^(X \E^ 



^(0)/^(0)-l 
^(0)/j^(0)-l 



(-1) 



T {x,t) = (X a \ 



^(0)^(0)- 



(o) 



(-1) 

Xo). 



|A ), 
|Ao>, 



(3.25) 

(3.26) 
(3.27) 
(3.28) 



Notice that in order to get the above relationships we have used the commutation rules for the corre- 
sponding sl(3) affine Kac-Moody algebra elements (B.4)-(B.ll) and (B.25)-(B.31), as well as their properties 
acting on the highest weight state |A ) (B.15). 

According to the solitonic specialization in the context of the DT method the soliton solutions are 
determined by choosing suitable constant group elements h in the eqs. (3.26)-(3.28). In order to obtain 
N— soliton solutions the general prescription is to parameterize the orbit of the vacuum as a product of 
exponentials of eigenvectors of the operators e; (ei = A vac , £2 = B vac ) defined in (3.6)-(3.7); i.e h — IlfL 1 e Fi , 
where [e; , Fi] = A'F^, , such that (i*i) m 7^ only for m < m,i, nii being some positive integer. The 
relationships between DT, solitonic specialization and the Hirota method have been presented in [15] for any 
hierarchy of integrable models possessing a zero curvature representation in terms of an affine KacMoody 
algebra. The DT method provides a relationship between the fields of the model and the relevant tau 
functions, and it explains the truncation of the Hirota expansion. The Hirota method is a recursive method 
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which can be implemented through a computer program for algebraic manipulation like MATHEMATICA. 
On the other hand, the DT method requires the computation of matrix elements as in eqs. (3.26)-(3.28) in the 
vertex operator representations of the affine KacMoody algebra. Actually, these matrix element calculations 
are very tedious in the case of higher soliton solutions. 

3.1 Free field boundary conditions and dressing group 

As a generalization of the constant NVBC (3.3) consider the free field NVBC 

lim tt± -+fif{x,t) = pfe a t x+b t t (3.29) 

x— y — CO J J J 

where pf , af, are some constants. Notice that (3.29) is a free field solution of (2.6)-(2.7) such that 6* = 
±(af) 2 T2Aj, for Aj = E,=i Pi Pi ~ \ = const. However, a direct application of the DT approach of 

[15] is not possible since the relevant connections A vac , B vac do not belong to an abelian subalgebra (up to 
the central term). So, in order to consider the NVBC (3.29) we resort to the dressing group [29] composition 
law of two successive DTs. In fact, the relevant connections can be related to the A vac ,B vac in (3.6)-(3.7) 
through certain gauge transformations generated by 0± such that 

vr c = e^z.ijv^e^z.t)]- 1 + d xi e 9 ± (x,t)[e 9 ± (x,t)}- 1 ; Vi = A vac , v 2 = B vac , (3.30) 

where i = 1, 2, a?i = x, x 2 =t,Vi = A vac , V 2 = B vac , and 9| G SL{3). One must have 

*(o) ^ £(0) = Q9+ ( X)t )^(0) = Q 9_ ( X)t )*(0) ffj v{ vac) = d xi ¥V[¥%\ (3.31) 

g is an arbitrary constant group element, and 

q9 = gXOgE^oX^ e s = e -xo e £ n>0 x-„ ; [D, Xn ]=nXn- (3.32) 

From (3.30) one has e - Xo E^e Xo = E®, so 

Xo = x^<-&+^^ 1+ , 2 +Exo^ 0) +xC (3-33) 

a=l 

2 2 2 

x-i = -^[/+(x,o-p+]^ 1) +E[/r(^o-^]^+E^-i^" 1) - ( 3 - 34 ) 

i— 1 i— 1 a— 1 

Therefore, a modified DT procedure can be implemented with [29] 

[e^]- 1 e^ = *(°)/i / [*W]- 1 = e»( a; ,i)* (o) M* (o) ]" 1 [0+(^i)]" 1 ; & = (3.35) 

through 

A' = Q h l 9 A^[Q h l 9 ]- 1 +d x e h ± 9 [e h ± 9 }- 1 ; B' = Q^B&^Q^ 9 ]- 1 +d t Q± a [Q± a ]- 1 . (3.36) 

So, the equation (3.35) can be used instead of (3.14) in order to derive general dark-dark solitons with free 
field NVBC. 
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4 Bright solitons and vanishing boundary conditions 

For simplicity we first apply the DT method to the VBC (3.2) and show the existence of bright-bright soliton 
solutions [30]. The connections A vac = E^\ B vac = E^ are related to the group element <J>o as 

A {vac) = x tfo[*o] _1 ; B {vac) =d t ^o[^o}~ 1 ; *o = e xE(1)+tE<2) . (4.1) 

In order to obtain soliton solutions the simultaneous adjoint eigenstates of the elements E^\ E^ play 
a central role. In the case at hand one has the eigenstates Fj and Gj in (C.l) such that 



= —<pj(x,t)Fj; ipj(x,t) = Vj (x + Vjt) (4.2) 



xE w + tE {2) , Gj = r)j (x, t)Gj , r)j (x, t) = pj (x + pjt) (4.3) 

4.1 j '^-component one bright-soliton solution 

Consider the product 

h = e a ^ F ^e b ^ G ™, (4.4) 

where j\ and 32 are some indexes chosen from {1,2}. Using the nilpotency properties of Fj 1 and Gj 2 (see 
Appendix C) one gets 

[^o^o 1 ] = (l + e-^a jl F jl )(l + e^b h G j3 ) (4.5) 

= 1 + e" a h Fj 1 + e"» b h G n + a n b h e"^ e"* F h G h , (4.6) 

with (pj 1 and r]j 2 given in (4.2)-(4.3). The corresponding tau functions become 

Tn - 1 + n-h-C- ■ p"Wip'»2 (7- - - Vjl Pn A. . (a 7) 

Wi - Ph) 

t+ = 8ij 2 b ja p ja e^ , rr = S itjl a n v h e"^ , (4.8) 

where the matrix element Cj 1 j 2 has been presented in (C.2). In order to construct one-soliton solutions we 
must have j\ = 32 = j in (4.7). Therefore one gets 1 

= hMH = a ^' e ' y ' i = j . (49) 

1 + aibiCue-^ev*' 1 1 + ai h Cu e"** ' - h K ' ' 

*± = 0, it j (4.10) 

Impose the relationships p* = —Vj, b* = —pScij isj, aj G C; so from (4.7) one has Cjj = — (j^) 2 - The 

equations (4.9)-(4.10) with the complexification (2.10) provide a solution of the CNLS system (2.11) 

ipi(x,t)=\ V^K^F (4.11) 
i + 3 



1 If ji 7^ ]2 in (4.7)-(4.8) one still has certain trivial solutions, since in this case Cj lt j 2 = implying tq = 1 
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where e x ° — l ff}V^*\^ , <Pj = Vj {x — ivjt) = (fjR + itpji- It must be 8 = +1, and the solution (4.11) possesses 
2 complex (aj, Vj) parameters plus the real coupling p > 0. 

The solution (4.11) for the j'th component is known as a 'bright soliton' in the context of the scalar 
non-linear Schrodinger equation (NLS). 

4.2 1-bright-bright soliton solution 

The main observation in the last construction of the j" 1 — component one-soliton is that it has been excited by 
the group element h in (4.4) such that ji = j 2 = j- So, in order to excite the two components of tyf (i = 1,2) 
and reproduce a bright soliton for each component let us consider the group element 

h = e aiFl e blGl e a2F2 e b2G2 , (4.12) 

where the exponential factors contain F's and G's of type (C.l). The tau functions become 



+ = bj pj e v * , f]j = pj (x + Pj t) + p 0j , j = l,2. (4.13) 



7 

tJ = cijVje^^ 3 , ipj = Vj (x + Vjt) + u j, j = l,2. (4-14) 
t = l + a^C^e-^e^ +a 2 b 2 C 22 e-^ 2 e^ 2 , (4.15) 

where Cjj = ( p P ^.yi ■ Let us set b* = —pSjCij, p* = —Vj = — fi,then 77* = — ipj = —ipi = —(^pm + 
i(fn), {vj,dj e C) in the relations (4.13)-(4.15). Therefore, the vector soliton solution of the 2-CNLS eq. 
(2.11) arises 

'^)J * ) VlR ~±U* + %)*», (4.16) 
Mx,t) J \A 2 J \ 2 J 

where A { = , . ^' ; e x « = p 12 / ■ 

This solution is valid for Y]j 8j\a,jVj\ 2 > 0. Notice that this 1-bright-bright soliton possesses, apart from 
the real p > 0, four arbitrary complex parameters, namely, 01,02, v\, v 2 . This solution in the case of mixed 
nonlincarity 5i = —5 2 = 1, may have a singular behavior if the sum \jY^=i $j \ a j v j\ 2 m l ne denominator 
of the expression of Ai vanishes, such that the soliton amplitude in (4.16) diverges. The N— bright-bright 
soliton requires the generalization of the group element in (4.12) as h — e a i F i e fc i G i... e ajvFlv e blvGjv . 



5 AKNS r (r = 1, 2), NVBC and dark solitons 

In this section we tackle the problem of finding dark soliton type solutions of the system (2.6)-(2.7). The 
associated 2— CNLS model (2.11) with nonvanishing boundary conditions has been considered in the frame- 
work of direct methods, such as the Hirota tau function approach (see e.g. [10, 31, 22, 32]), and recently 
in the inverse scattering transform approach [12, 13]. In the last approach the relevant Lax operators have 
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remarkable differences and rather involved spectral properties as compared to their counterparts with van- 
ishing boundary conditions (see [34] and references therein), e.g. in the NVBC case the spectral parameter 
requires the construction of certain Riemann sheets [34, 35, 36, 37]. So, it would be interesting to give the full 
Lie algebraic construction of the tau functions and soliton solutions for the system (2.6)-(2.7) with NVBC. 



5.1 AKNSi: N-dark solitons 

For simplicity, firstly we describe the entire process for the system (2.6)-(2.7) with just two fields i $ ± . So, 
let us consider the sl(2) affine Kac-Moody algebra in the Weyl-Cartan (WC) basis 2 



[H (m) ; H { n)] = ™ Sm+nt0C , \H^ m ) , 4 n) J = ±4 m+n) , [E { ; n) , = 2H^ + m5 m+nfi C. (5.1) 

In order to study a NVBC for the system sl(2)-AKNS consider the Lax pair 

A = + *+4 0) + *~£i 0) + 0i C, (5.2) 

B = E {2) + V+E^ + + 

d^+E^ - d^-E^ - 2 (*+*- - p+p-) H {a) + <t> 2 C, (5.3) 

where arc the fields of the model (p ± = constant). In this case one considers the generator 

= H~( l \ The 4>i and 4>2 are introduced as auxiliary fields. Therefore the equations of motion suitable to 
treat NVBC become 

= fl£tf+- 2 (*+*-- p+p') tf+ (5.4) 
9 t *- = -a 2 *" + 2 -p+p-)^-, (5.5) 

&t<h - d x <}>2 = 0. (5.6) 

The AKNSi model (5.4)-(5.5) is recovered from the AKNS 2 extension (2.6)-(2.7) simply by setting to 
zero the additional fields. In [33] it has been introduced a complexified version of the system (5.4)-(5.5) 
(for P+ = 0) and presented its reduction to the focusing and defocusing NLS system, as well as its soliton 
solutions. In fact, the sZ(2)-AKNS model (5.4)-(5.5) as well, through a particular reduction, contains as 
sub-model the scalar defocusing NLS system 

id t ^ + dl^-2 (|Vf-p 2 )^ = 0. (5.7) 

This equation is suitable for treating nonvanishing boundary conditions (NVBC) [37, 34, 35] 

{p, x — > — oo . . .„ 

; e = e' e / 2 , pel. (5.8) 
pe 2 , x — » +oo 

In fact, this boundary condition is manifestly compatible with the eq. of motion (5.7). The defocusing 
NLS equation (5.7) with the boundary condition (5.8) is exactly integrable by the inverse scattering technique 



2 This basis differs from the Chevalley (Ch) basis in (A.1)-(A.4) by the rescaling of the generator — ► iH^^. 
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[38] . This model has soliton solutions in the form of localized "dark" pulses created on a constant or stable 
continuous wave background solution. 

In the system of eqs. (5.4)-(5.5) consider the transformation 

t — > —it, x — > — ix, (5.9) 
* ± -> i* ± e =F2 , (5.10) 

where the factor e T2 is introduced for later convenience. Furthermore, the identification 

ip = *+ = (*-)*, (5.11) 

where the star stands for complex conjugation, such that p+p - — > — p 2 allows one to reproduce the defocusing 
NLS equation (5.7). 

Let us take as vacuum solution of (5.4)-(5.5) the constant background configuration 

* ± = p ± = pe T2 , </>i,2 = 0, , p, e = const. (5-12) 

Therefore the gauge connections (5.2)-(5.3) for the vacuum solution (5.12) become 

A vac = £i = E {1) +p + E i ° ) +p-E W (5.13) 
B vac = e 2 = + p+ E { + ] + p- E W (5.14) 

We follow eq. (3.8) in order to write the connections in the form 
with the group element 



,(°) _ „xsi+te 2 



= e 1 

^ nvbc — c 



exp 



x(EW + p + E f + p-E^) + t(E^ + p+E ( ^ + P -E { H (5.16) 



Let us emphasize that this group element differs fundamentally from the one associated to the case with 
VBC. First, the difference originates from the constant boundary condition terms added to the relevant 
vacuum gauge connections. However, they must be related since lim *^ bc — > ^o, where *o is the group 
element in the VBC case (4.1). Second, the vacuum connections £1,2 in the DT procedure will require another 
eigenvectors under the adjoint actions in analogy to eqs. (4.2)-(4.3), as we will see below. 

Another proposal for this group element *^ bc has been introduced in [39]. There it has been considered 
a group element inspired in the inverse scattering method, possessing double- valued functions of the spectral 
parameter A. This fact motivates the introduction of an affine parameter to avoid constructing Riemann 
sheets. This construction involves some complications when used in the context of the DT method. The 
main difficulty arises in the computation of the matrix elements associated to the highest weight states in 
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order to find the relevant tau functions. However, the group element given in (5.16) will turn out to be more 
suitable in the DT procedure. In fact, a similar group element has been proposed in [40] to tackle a NVBC 
soliton solutions in the so-called negative even grade mKdV hierarchy. 

The DT procedure in this case follows all the way verbatim as in section 3 adapted to the si (2) case. It fol- 
lows by relating the relevant connections (5.2)-(5.3) with the connections in eqs. (5.13)-(5.14) corresponding 
to the vacuum solution (5.12). 

Next, we will write the fields in terms of the relevant tau functions. According to the development 
in section 3 these functions will be written as certain matrix elements in an integrable highest weight 
representation of the affine Lie algebra sl(2). So, it is a straightforward task to get the following relationships 



*+ = p + H and *~ 

to 



_ T 

p 

TO 



where the tau functions , r are given by 

T+ EE (A |£« 
T~ EE (X \E^ 
To EE (A | 

nvbc § iven b y ( 5 - 16 ) 



>(0) 
■ nvbc 1 



(0) 

nvbc' 



>(0)-l 

" nvbc 



(0)-l 




nvbc 


(-1) 


(o)-r 




nvbc 


(-1) 


|A ), 


J(o) 





|Ao) 
|Ao) 



(5.17) 

(5.18) 
(5.19) 
(5.20) 



with the group element 



5.1.1 AKNSi and 1-dark soliton of defocusing scalar NLS 

Let us consider the group element 

h Q =e a«v«( 7 ,p ± ) } ( g = l,2) (5.21) 

where a q is a constant. The vertex operator V q satisfies V q — 2V q (see below) where V q is defined in (A. 10). 

Next, one must look for the eigenvalues of s\ and e 2 , respectively, in their adjoint actions on the vertex 
operator V q . In order to achieve this, let us notice that these Weyl-Cartan basis elements £1,2 are related to 
the corresponding Chevalley basis elements i\ and e 2 in (A. 12) through i n — 2 e n , (n = 1, 2), provided that 
p ± = 2p ± and the relationship between the Weyl-Cartan and Chevalley basis H^ c — is considered 

for the sl(2) element H^ n \ So, the vertex operators are related by V q = 2V q . Then from (A. 11) one has 

[ Sl ,V q }= 7 V q , [e 2 ,V q ] = (-l) q - 1 1 ( 1 2 -4p + p-) 1/2 V q , g = l,2. (5.22) 

Notice that these eigenstates V q exhibit a more complex structure in comparison to their counterparts 
Fj and Gj corresponding to the VBC as in (4.2)-(4.3). Therefore one has 

[xe l +te 2 ,V q { ll p ± )} = [ 1 (x + (-l) q - 1 vt)]V q (^p ± ), (5.23) 

where v = \J~f 2 — 4p+p~. Denoting ip q = 7 (x + (— l) q ~ 1 vt) one can write 

= ex-pie^aqV 1 ) (5.24) 
= l + e v "a q V g , (5.25) 



* (0) h^- 1 

nvbc nbvc 
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Vector Soliton^, ¥") 
4h 




Figure 1: The real vector (^ + , "J ) 1-soliton profile of the AKNSi model for q = 1 
in (5.26). They were plotted for 71 = 8, p+ = p~ = 3.9, logcW = 8. The solid and 
dashed lines correspond to \f r+ and V? - , respectively. 



where we have used (V q ) n = 0, for n > 2, that is, V q is nilpotent (A.17)-(A.18). 

So, in order to find the explicit tau functions in (5.18)-(5.20) it remains to compute the relevant matrix 
elements. Using the properties presented in the Appendix A and the eqs. (3.26)-(3.28) one gets the tau 
functions r° q) = 1 + e v «; rf q) = s ± ^p ± e^", where the matrix elements = (A Q | V q |A ) , s ± («> = 
{p ± )~ 1 (X \E^V q |A ) have been considered. They are given in the Appendix A, eqs. (A.13)-(A.14), 
respectively. So, one has S 2 J^ — 4p + p - ( 7 i ( — l) 9_1 \/7 2 — ■ Without loss of generality we can 

assume 7 > 0, so, the matrix element (A | V q |A Q ) from (A. 13) allows us to determine the sign of c^ q \ which 
is completely fixed by sign[c^ 9 ^] = e q sign[a^ qS> ]. We will see below that the signs of the free real parameter 
c*- 9 ' determine two different type of solutions. 

Consider > 0. This case follows if as> q ' and e q have equal sign. Then, the relations (5.17) supplied 
with the above tau functions provide 

^\x,t) = p ± [l- w ] - P -^ r tanh{ 7 [z + (-iy- 1 vt]/2+ -L— }. (5.26) 

These solutions are the vector 1-soliton of the AKNSi model. This solution possesses 4 real parameters, i.e. 
p ± , 7, . Notice that the value of the index q determines the direction of propagation of the soliton. The 
Fig. 1 displays the component solitons of this vector 1-soliton. The true 1-dark soliton of the AKNSi model 
is obtained as the product 

[*- (9) * + (q) }(x, t) = p+p- - ^ S ech 2 { 7 [x + (-l) q - 1 vt}/2 + (5.27) 

It has been plotted in Fig. 2 for the case q = l. 

Next, consider c( q ) < 0. This case follows if a^ q ' and e q have opposite sign. Similarly, the relations (5.17) 
and the above tau functions provide 

s ±(<3) o ± s ± ^ , loa\c ( - q '>\ 

^(x,t) = P ± [1- - P -^ r coth{ 7 [x+ (-!)«- V]/2+ -^—4- (5.28) 
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-t 0.0 

Figure 2: The 1-dark soliton profile of the AKNSi model emerges as l I' _, J ,+ . This 
is the soliton traveling to the left, the case q = 1 in (5.27). It was plotted for 
7l = 8, p+ = p~ = 3.9, logcW = 8. 

These solutions are the singular solitons of the AKNSi model and, to our knowledge, they have never been 
reported yet. As above the value of the index q determines the direction of propagation of the soliton. They 
possess 4 real parameters, namely, p^, 7, c^ q \ Similarly, let us consider the product of the components 

(«)] (x,t) =p+p- - (-l)^ 1 ^csch 2 { 7 [x + (-l) q - 1 vt}/2 + l ° 9 1 } (5.29) 

These functions are displayed in Fig. 3. Notice that these solitons are singular for [7 [2; + (— l) 9 ~ 1 vt]/2 + 

Some comments are in order here. 

First, the 1-dark soliton solutions (5.27) and the singular solitons (5.29), respectively, have the same 
boundary values in x — > ±00, i.e lim 'I' - ^ q \x, t)\P + ^ (x, t) — > p + p~- In the both types of solitons 

I a; I — ^+oo 

this behavior is in contrast with the different boundary values assumed by the corresponding components 
in the both limits. For the 1-dark soliton components one has lim , J ,± — > p ± , whereas lim — > 

X— > — 00 ' x— >+oo 

_i_ 7 , 

p [1 — (7 ± v7 2 — 4p+p~)l. This fact is shown in Fig. 1. 

2p+p 

Second, the center of the 1-dark solitons are given for {j[x + {-l) q - 1 vt]/2 + = 0. In fact, in this 

point the functions ^^f + ^ q '\ C enter = P + P~ — 7 2 /4 (for q = 1,2) take the smallest intensity. So, their 
intensity dips are controlled by the value of the parameter 7. 

Third, the solution (5.26) in the case q = 1 allows the imposition of the conditions (5.11) in order to 
satisfy the defocusing NLS equation (5.7). So, let us make the transformation (5.9) and the complexification 
71 — > 17, p — > ip, v 1 — > v 7 — y/ Ap 2 — 7 2 ; (2p > 7). Choose a± = —ij/v. Moreover, considering the 
transformation (5.10) one gets the functions satisfying the complexification condition (5.11). Therefore, the 
function Tp(x,t) = —ie 2 ^/ + , where ty + = ie~ 2 {p + ( -f- 3 ^ ) y0m^jt } 1 vi x i~t) — l( x ~ %t) + 7o, obtained in 
this way is a solution of the defocusing NLS equation (5.7). In fact, the function tp can be written as 

ip(x,t) = e l(Q+7r/2) {^tanh[7(x - v„ ( t)/2 + 7o /2] + 1- ip e - ia }, a = tan" 1 ^), (5.30) 
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Figure 3: The singular solitons of the AKNSi model. The solid lines travel to 
the left (q = 1) and the dashed lines to the right (q = 2). They were plotted for 
t = -0.07, 7 = 10, p+ = p- = 2, logical = 10. 

which is the dark soliton solution of the defocusing NTS equation (5.7) (see e.g. [31]). Notice that this 
solution satisfy the NVBC lim \ip(x,t)\ 2 ^ p 2 . 

x— >±oo 

Fourth, similarly, one can consider the complexification condition (5.11) in the case q = 2. In this case 



one gets i & ± = ie^ 2 f 



T ^e ±i,5 tanh[(7a; + 7« 7 t)/2] 



2p 



, S = tan 1 {v 1 /^j). However, they do not 



satisfy the condition (5.11), hence these functions would not be solutions of (5.7). 

From the previous constructions it is clear that the N-solitons of the types dark, singular or mixed dark- 
singular can be constructed by choosing convenient signs of the parameters a? and e q associated to the 
group element 



1,2 



(5.31) 



5.2 AKNS 2 : N-dark-dark solitons 

Consider the non-vanishing boundary condition (NVBC) for the system (2.6)-(2.7) in the form (3.1) and 
its associated constant vacuum solution (3.3). The vacuum connections (3.6)-(3.7) associated to this trivial 
solution commute [A vac , B vac ] = if one takes Ox = ^2 = |(Pi Pi + ptPz)- Consider the notation 

A vac = El , B vac = El. 

1-dark-dark soliton. Let us consider the eqs. (3.26)-(3.28) and choose the group element h q = e aqwq ( k ) ; 
a q and k q being some parameters, where the vertex operator W q (k) is defined in (C.7), such that 

[ei,W q ] = k q W q ; q=l,2; (5.32) 

[e 2 , W q ] = w q W q (5.33) 

where 

w q = e q k q ^{k q Y - Ap+p^ - 4p+pz, e , = (-1)9-1 (5 34) 
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So, one has 

[xei + te 2 , W q ] = (k q x + wH)W q (5.35) 

The vertex operator W q (see (C.7)) is associated to the 1-dark-dark soliton and it can be shown to be 
nilpotent (C.8). So, using (3.26)-(3.28) the following tau functions correspond to the element h q given above 



t, 



(a) 1 _1_ r {l) p k q x+w q t 







l + c^ q 'e k x+w \ (5.36) 
T f {q) = p f s ±(i) e k^+^t^ i = 12 . (537) 

where = a q < X \W q \X >, s ±(g) = {pf^cfl < A |£^. W q \X >. The above tau functions provide a 
1-dark-dark solution for the fields ^>f of the system (2.6)-(2.7), such that the following relationships hold 
between the parameters 

C W = ^ 8 ■ (fc*) 2 = -if s ' y^kPkPkl w , = s +s (fc 9)2. (5. 38) 

s+ (<?) - s~ («) s+w)s~w) s+ (9) — s~ w) 

We will concentrate on the regular solutions of the AKNS2 model since they will give rise to the dark solitons, 
as it was clear from the AKNSi construction. So, in what follows we will require > 0. Then, the relations 
(3.25) provide 



± (?) _ l J j 



(l + (y«) ±1 ) + (-l + (y q ) ±1 )tanh{(k q x + w q t + logc^)/2} , (5.39) 



where j = 1,2; y q = s s _ ( ^ . The condition > requires that y q > land y q < 0. This vector soliton 
possesses 6 real parameters, i.e. p^, k q , c^ q \ 

Taking into account (C.7) one can compute the relevant matrix elements, which define s ± ( q \ in order to 

get 

y q = V . = ■ (5.40) 

k q ~e q ,J(k<iy-4j: 3 p+p- 

Notice that for any value of the index q, if J2j Pj Pj > the last relationship implies y q < and if 

Ej pjpj < thcn v q > °- 

Some remarks are in order here. 

First, the solitons associated to the pair of components (vE'^" , \1/^") and (v&T/ , \E'^~), respectively, are 
proportional, so they are degenerate and reducible to the single dark soliton of the sl(2) AKNSi model (5.26). 
Whereas, the solitons associated to the pair of components (tyf , respectively, arc not proportional, 

so they arc non-degencrate presenting in general different degrees of 'darkness' in each component. In this 
way the solution in (5.39) presents a non-degenerate 1-dark-dark soliton in the components, say (\E f j l " , "J/^ ) 
of the AKNS2 model. A nondegenerate 1-dark-dark soliton in the defocusing 2-CNLS has been given in 
[22, 10] through the Hirota's method. The inverse scattering [12, 13] and Hirota's method [10] results on 
N-dark-dark solitons in the defocusing 2-CNLS model have presented only the degenerate case. 
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Second, as in the AKNSi solution (5.26), one can show that the 1-dark-dark solitons (5.39) allow the 
imposition of the conditions (2.10) for the pair of fields (tyf , ^~), i = 1,2, in order to satisfy the 2- 
CNLS equation (2.11). In fact, consider the parameters y q — z q f\ and a complexification procedure as 
in (2.10) provided with a convenient set of complex parameters in (5.39), k q — > ik q , pf — > ipf such that 
y q = e 2 ^", z q = ^^J^J2j Pj Pj ~ (k q ) 2 , 4>q being real, so the equation (2.10) can be satisfied if 

( P +y = pS jp -. (5.41) 

Notice that, as it was shown in the AKNSi case and its 1-dark soliton (5.30) particular reduction, in order 
to get 2-CNLS 1-dark solitons one must consider q = 1 in the both components (^f , Therefore, the 

relationship between the parameters pf determine clearly if these solutions will correspond to the defocusing 
(Sj = —1) or to the mixed nonlinearity (5i = —S 2 = ±1) 2-CNLS system, respectively. 

Third, let us discuss the extension of our results to the case of free field NVBC (3.29). In order to compute 
the relevant tau functions analog to the ones in (3.26)-(3.28) one must consider the expression (3.35) instead 
of (3.14). So, the tau functions become 

f+ 9 (x,t) = (XolE^e^x^^h^}- 1 ]^) e"* (5.42) 

= (pf(x,t)-pf)f 9 (x,t) + (X \ [e-x-E^l e^^hi^}- 1 ^) (5.43) 

rr 9 (x,t) = (X \E { ^e 9 _(x 7 t)^ a '>h[^}- 1 \\ ) e"X, (5.44) 

= (-pr(x,t)+pr)ri(x,t) + (A c | [e^E^e**] ^h[^]^ \\ ) (5.45) 

f°(x,t) = (\ \ * W/i^W]" 1 |A D ) . , (5.46) 

where the group element 9^ (x, t) is defined in (3.31)-(3.32) and x{x, t) is an ordinary function. We have used 
Q 9 + from eq. (3.32) and the decomposition (3.33) for xo> as well as the properties (B.15)-(B.16) of the highest 
weight representation. The tau functions (5.42)-(5.46) exhibit the modifications to be done on the previous 
eqs. (3.26)-(3.28) in order to satisfy the NVBC (3.29). These amount to introduce the group clement 6^ 

- + g - — g 

and the factor e~ x . Then from (5.42)-(5.46) and the analog to (3.25) *+ = pf + and *r = pT - T±^- 
one can get 

*+ = p+(x,t) + ^; *r=pr( x ,t)-?4- (5.47) 

TO To 



ft 



(A | \e-*°E { \e*°] ^h^]- 1 \\ Q ) ; f 9 = f = (A | ^h^]- 1 |A > (5.48) 



In fact, the relations (5.47) can be formally obtained from the ones in (3.25) by making the changes 
pf — > pf(x,t), rf — > ff, tq — > fo. The soliton type solutions can be obtained following similar steps 
as the previous ones; i.e. for 1-soliton choose h q = e aqwq ( k ) and the adjoint eigenvectors of £1,2 become 
precisely the vertex operator W q as in (5.32)-(5.33), then the expressions (k q x + w q t) of the tau functions 
(5.36)-(5.37) do not change. 

Next, let us describe the general properties of the solutions (5.47)-(5.48). i) Notice that the form of f 
remains the same as the one in (5.36). ii) The modified DT and the related tau functions f± (5.48) will 
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introduce some new parameters into the tau functions t~~ in (5.37). In (5.37) the factor pf must be replaced 
by pf(x,t) and s ± ^ in general will depend on the index l v and the space-time, i.e. s ± ^ — >■ s ± ( q '^ (x,t). 



W q \X >. In the case of 



The explicit form one gets from s ± ^ q ' l \x,t) — (pf [x, t)) _1 a 9 < A | 
1-dark-dark-soliton it is required an element \ such that s ± (x,t) = s ± ^ xf s ±( - q \ xf ^ being 
constant parameters. So, under the above modified DT one has the general form of the 1-dark-dark soliton 



± (q) _ P] e 



±„a :t x+b :t t 



1 + (y«) ±1 ) + (-1 + (yf)* 1 ) tanh{(k*x + wH + logcf)/2} , (5.49) 



i 2 
where j = 1, 2; = p-^J, and 

2 1 

&± = ±{aff T 2£p>T - -CM)]; M ^ /3 2 .0. (5.50) 

3=1 

The solutions (5.49) are the general 1-dark-dark soliton components of the AKNS2 model. The index 
q = 1,2 refers to the solitons traveling to the left and right, respectively, and this solution has 13 real 
parameters, namely, pf, k q , c^\af, fii, S^2 ; i.e. additional 7 real parameters as compared with the solution 
in (5.39). A complexified version of this soliton for q = 1, which is a solution of the 2— CNLS model, has 
recently been reported in [4]. Moreover, the general form has also been reported in [14] for AKNS r using 
another approach. As in the AKNSi case (5.27), it is clear that the 1-dark soliton profiles can be recovered 
by plotting the functions [*T (<?) *+ (q) ](x, t), (j = 1,2; q = 1,2). 

Moreover, similar steps can be followed in order to get the singular solutions of the AKNS2 model related 
to Cj 1 ^ < 0. They will give rise to singular solitons as in the case of the AKNSi construction (5.28)-(5.29). 
It is also possible to construct mixed dark-singular solitons by conveniently chosen the positive and negative 
values of the parameters of each component of the vector (^f)- 

2-dark-dark solitons. In order to get the two-dark-dark soliton (N = 2) solution one must choose the 
group element h\ = e a ^ wq ( k ^) e a ^ wq ( k ^) _ Then following similar steps one has 

r = l + Cl e klX+Wlt +c 2 e k2X+W2t + Cl C2d e klX+Wlt e k2X+W2t , (5.51) 
rf = pflsf e klX+Wlt +s±e k2X+W2t + sfs±d ± e klX+Wlt e k2X+W2t ], i = 1,2 (5.52) 

where 

2 




d = | V % lS2 ;d ± = - ^\. 44 \ d ; (5.53) 

' {s 1 -s+)(s 2 -sj) 



S t S 2 



i^- ; (fc„) 2 = - ( ' S " Sn )2 | Efc=1 Pk Pk ] ; w n = k 2 n , (n = l,2) (5.54) 

Notice that the above tau functions related to two-dark-dark soliton must require Cj > 0. This two-soliton 
has 8 real parameters, i.e. p- , fci, k 2 , c\, c 2 . The above process can be extended for N— dark-dark soliton so- 
lutions, in that case the group element in (3.26)-(3.28) must take the form h\ = e a ^ w i e a2W % {k2 \..e aN w N( k ") . 
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The relevant tau functions follow 

TO = 



; m,n = 1,2,3. ..N (5.55) 

iVxJV 



pf ( <y mn + + cmn e^ 3 ^' - ll , i = 1,2. (5.56) 

^ NxN J 



where | |jvxiv stands for determinant, c mn = ,„- m " *g " + . ( ; = = ^7== ) , and w„, fc n are the same 



2 



as in (5.54). Notice that c mn — c nm , c nn — 0. We have omitted the index (q) in all the parameters above. 
The N-dark-dark soliton (5.55) possesses 2A+4 real parameters, i.e. Pj , k m ,c m . Let us mention that we have 
verified these solutions up to N — 3 using the MATHEMATICA program. The above solutions associated to 
the two (5.51)-(5.52) or higher order (5.55)-(5.56) dark-dark-solitons are in general non-degenerate. Notice 
that in the CNLS case the two— and higher— dark-dark solitons derived in [10] are actually degenerate and 
reducible to scalar NLS dark solitons. So, regarding this property our solutions resemble to the ones recently 
obtained in [4] for the r— CNLS system and in [14] for the AKNS r model. 

Let us remark that the free field NVBC (3.29) requires the introduction of more parameters into the 
N-soliton tau functions (5.55)-(5.56) through the modified DT. So, following similar steps to get (5.49) one 
has: i) the parameters entering the tau function To in (5.55) remain the same, ii) in (5.56) the factor 

must be changed to pfe aj x+hj *• hi) the parameters entering the tau functions in (5.56) must be 
changed as — > s^xf- This general N— soliton will have, in addition to the Oi ; 2 parameters, 3A^ + 8 real 
parameters, i.e. p^ , k m , c m j, aj . 

Similar constructions can be performed to get the singular and the mixed dark-singular N-solitons. In [4] 
the general non-degenerate N-dark-dark solitons in the r-CNLS model with defocusing and mixed nonlinearity 
have been reported in the context of the KP-hierarchy reduction approach, and in [14] the bright and dark 
multi-soliton solutions of the AKNS r system have been addressed in the algebro-geometric approach. 

5.3 Dark-dark soliton bound states 

So far, reports on multi-dark-dark soliton bound states in integrable systems, to our knowledge, are very 
limited. Recently, it has been shown that in the mixed-nonlinearity case of the 2-CNLS system, two dark-dark 
solitons can form a stationary bound state [4]. Then, in order to have multi— dark— dark soliton bound states 
in the N-soliton solution (5.55)-(5.56) the constituent solitons should have the same velocity, i.e. denoting 
y n = j^, (y n ^ [0,1]) then ^ = y y n ^_\ k n = v (assume v > 0) for certain soliton parameters labeled by 
n = 1, 2, .... We show that the signs of the sum pj determine the existence of these bound states, for 

the positive sign bound states can be formed, whereas for the negative sign they do not exist. 

First, consider Y^iPtPi > 0> tncn f rom (5.54) it follows that k n — j (X^ Pi Pi ) 1 ^ 2 f° r 2/n < . 
Therefore, in order for multi— dark— dark soliton bound states to exist the equation 



= c>0 (5.57) 



21 



must give at least two distinct positive solutions for \y n \. In fact, one has the following two solutions 



c 2 yi6+? yi 2c 2 + c 4 + A 
|W, 2 | = 1 + T + — — ± W2 > (5-58) 

where A = 64c ^f^ c ^ c ■ These exaust all the possible solutions with the condition (5.57). The right hand 
side of (5.58) with either + or — signs provides 1 2/1,2 1 > for any c > 0. The other possible solutions 2/3,4 do 
not satisfy the above requirements, since they possses a term of the form \J\2c 2 + c 4 — A, which is imaginary 
for any positive value c. Therefore, two— dark— dark— soliton bound states exist in the sl(3)— AKNS system, 
and three— and higher— dark— dark— soliton bound states can not exist. These results hold for any value of 
the index q. Notice that in order to reduce to the 2-CNLS system the parameter relationship (5.41) must 
be satisfied. So, from (5.41) one has J2iP^P7 = 77 1^77 — ^ T7~] < ^ ( remcm bcr that under pf — > ipf 
the sum J2i Pt PT reverses sign); so, it is possible if Si = —5 2 (e.g. Si = —S 2 = —1 for \p^\ > \p 2 \), this 
case corresponds to the 2-CNLS with mixed focusing and defocusing nonlinearities. Thus, the two-dark- 
dark soliton bound state solution we have obtained here corresponds to the Manakov model with mixed 
nonlinearity. The case Si = = 1,2) defines the defocusing Manakov model which does not support 

multi-dar-dark-soliton bound states [4] . 

Second, the condition pfp~ < will provide only a single positive solution for y q > in the equation 
^— = v > 0. So, one can not obtain two or more solitons with the same velocity and therefore bound states 
in this case are not possible. 



6 Mixed boundary conditions and dark-bright solitons 



One may ask about the mixed boundary conditions for the system (2.6)-(2.7), i.e. NVBC for one of the field 
components of the system, say ^f, and VBC for the other field component, 4^. So, we will deal with the 
mixed boundary condition (3.4). Moreover, taking into account sequence of conditions as in (2.10) 

[*+]* = -pJi*:- s-MVi, (6.1) 



t 



-i t. 



the system (2.6)-(2.7), for (/?j.p) being real, can be written as 

1 



yj k = 0, fc = l,2. 



(6.2) 



Precisely, the system (6.2) in the defocusing case ( Si = S 2 = —1) and possessing the first trivial solution 
(3.4), i.e. pi 7^ 0; p 2 — 0, has been considered in [26] in order to investigate dark- bright solitons which 
describe an inhomogeneous two-species Bose-Einstein condensate. The system (2.6)-(2.7) with the mixed 
trivial solution (3.4) can be written as 



Z, 

±d 2 x ^ T 2[J2*l*i-pipi}*i> 



2 



(6.3) 
(6.4) 
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From the previous sections we can make the following observation: the vacuum connections relevant to 
each type of solutions must exhibit the fact that dark solitons are closely related to NVBC, whereas bright 
solitons to the relevant VBC one. Since and f2i are some parameters satisfying 2f2i + ^2 = 2ptPi> f° r 
simplicity we will assume O2 = 0, Q\ = pfpi in (3.6)-(3.7). Therefore, the connections (3.6)-(3.7) for the 
mixed (constant-zero) boundary condition (3.4) take the form 

i- ee ^^W+pf^+ft^, (6.5) 
B vac ee ^^E^+p+E^+p-E^l (6.6) 

The relevant group element is given by 

VQe ■ (6-7) 

In order to construct the soliton solutions we must look for the common eigenstates of the adjoint action 
of the vacuum connections (6.5)-(6.6). So, one has 



i,r ±/Ja ] = ±fc ± r ±fe , (6.8) 



[4,r ±j8a ] = i^r^, w ± = (k^f-p+p- (6.9) 

where the vertex operators r± ( g 2 (fc ± , pf) are defined in (C.3). We expect that these vertex operators will 
be associated to the bright-dark soliton solutions of the model. 
Let us write the following expressions 

[se^ + tej^ , I> 2 ] - <p + (x,t)Tp 2 ; y + {x,t) = k+x + w+t; (6.10) 
[xe\ + te 2 ^ , T-fa] = -ip-(x,t)T-i3 2 ; f~(x,t) = k~x + w~t. (6.11) 

The 1-dark-bright soliton is constructed taking h in (3.26)-(3.28) as 

h = e 7+r +f2 e 7_r -f2 , 7 ± = constants, (6.12) 

where the vertex operators in (C.3) have been considered. Notice that due to the nilpotency property of the 
vertex operators, as presented in the appendix C, the exponential series must truncate. So, replacing the 
group element (6.12) in (3.26)-(3.28) one has the following tau functions 3 



Ti = 



7+7" <A |4^r +& r_ & |A > (6.13) 
ft = e^V < X O \E^T ±02 \X O > (6.14) 
f = l + e^-^7+7- <X \r+p 2 T-p 2 \X > (6.15) 



3 There exist other eigenstates [e^ , V|J = A V£ ; [e 2 ^ , VjjJ = (-l) 9 " 1 A(A 2 -p 2 ,) 1 / 2 Vg« , where p 2 = 4p+p~ and (q = 
1, 2) is given in (C.5). However, these eigenstates are related to purely dark-solitons for the first component tyf, e.g. if one takes 
the group clement h = e^Pi , it will not excite the second component since the matrix elements of type < \ \E^p^ |A > 
vanish. 
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These tau functions resemble the ones in (5.36)-(5.37) for the 1-dark soliton (5.39) and the eqs. (4.7)- 
(4.8) for the 1-bright soliton (4.9) components, respectively. The matrix elements in (6.13)-(6.15) can 
be computed, and will depend only on the parameters k^jpf. So, one has six independent parame- 
ters 7 ± , fc* , p\ associated to the these tau functions . Let us write the tau functions in the form fo = 

1 + Coe k + x+™ + t e -k-x-w-t^ f ± = a ± e k+x+w+t e -k-x-w-t^ f ± = b ± e ±k ± x±w ± t^ where w ± = ( fc ±)2 _ 

PtPT\ c o = a+p °- + _V p + ; £ = ^;b+b- = - l)(k+k- +p+p-). So, using these tau functions in 

(3.25) one has 



a 



± 



'2co' 



(fc+ — k )x + (w + —w)t 



tff = pf±( — ) 1 + tanh^ '—- -i '- (6.16) 



2 



*± = -4= e ^[(fe + + fe-)^(^ + ^-)*l SC ch (fc+ = fc ' )x ± (W+ = W ' }t . (6.17) 
2vco 2 

This is the 1-dark-bright soliton of the si (3) AKNS model (6.3)-(6.4) (for tt 2 = 0). Notice that this solution 
has six independent real parameters, say a~ , b~, k ± , p l . 

The construction of the 2-dark-bright solitons follows similar steps. The group element 

h = e 7i + r +f ) 2 (fc+) e 7rr_ / 3 2 (fe7) e 7+r +fe (fc+) e i2 T -n 2 { - k 2) (6.18) 
does the job. We record the relevant tau functions 



T = 1 + 

e k+x+w+t e -k-x- w -t _|_ Co gfef x+tuftg-fcrx— uj^tgfejx+wjtg-fejx— w^t (6.19) 



2 



m,n=l 



-± _ ^ ^ a ± e *;+a:+u)+ *g- fe „ t _|_ ^igfc+x+uj/" *g- fe i x— ti^ tgfcjx+u)^ i g -fc 2 x--tu 2 t (6 20) 
m.n—1 

2 2 

f t = ^b^e^^ + e ±k ^ x±w ^ t e ±k ^ x±w ^ t [J2 ( l^ eTklxTwl % (6-21) 



m— 1 



l ± /; ±\2 + — a i,„ a ™„ a J,„ k i,pt dfd, df kfkfpf 

where w= = (/c*)" 2 - p{ p l ; c m „ = + "L" m _" + ; = c = ,+ - 1 ' + ; -jt = ,i, 2 - - ; 

Omnd — a mn pj a m „ fen Pi "i Pi — «i Pi d 1 fej fe 2 

°12( fc ] f -' c 2")( fe l hfc 2"+P] F Pr) I a22( fc 2~- fe 2 )( fe 2 fc 2 +Pl~/» r ) _ Q. ?,+ /,- — ( k t ~ k i )( k t k i + Pl Pl ) „~ . 

a" (fe+ - k~){ktk- +p+p-))J aji (fe " - fe+ ) (fej fe+ +p+ p" ) ' ^ ^ "" ' _ fc-p"_ " " «' 

,_ _ a~a 2 ~ 1 (fc+-fc+) 2 (fc~-fc 2 ~) 2 (fc+fc+-fc~fc 2 ~)(fc+fc++p+p~)(fc~fc 2 ~+p+p 1 ~) i 6~ _ a'fcj- (fc+ -fc~ ) (fc+fc^ +pj>£ ) . 

i " (fe+-fcr)2(fcr-fe+)(fc+-fc 2 -)(fc+-fc 2 -) 2 pr(fc+fer+pfpr)(fcJfc2"+p^pr)' &r a-fc 2 - (fcf-fcr)(fcffcr+pfcr)' 

«f = T f lfJ ( P f) 2 (kt-wf-kti) ^^ + PtP^ a Unf> m t =nr=3-i; mT = n+ = i. 

Notice that this solution possesses 10 independent real parameters, namely, a^ nm , a>i 2 , b± , fc^, pf, (m — 
1,2). These tau functions resemble to the ones in [24] provided for the 2-CNLS model. The generalization 
to N-dark-bright solitons requires the group element h to be 

h = e fi f r+^ 2 (fe+) e 7 1 "r_ / 3 2 (fe-) e 7+r + 3 2 (fe+) ^2^^^^-) ^^ e j+r +l32 (k+) e7 ^r_^ 2 (fe^)^ (6.22) 
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7 Generalization to AKNS r (r > 3) model 

The procedures presented so far can directly be extended to the AKNS r (r > 3) model for the affine Kac- 
Moody algebra sl(n) furnished with the homogeneous gradation. According to the construction in [20], in 
this case the equations of motion will describe the dynamics of the fields ^ (j = 1,2, ...,r; r = n — 1) 
associated to the generators E±p , where the /3j are the positive roots defined by fij = aj + Oj+i.... + a r 
(ctj =simple roots). The outcome will be the eqs. in (2.6)-(2.7) with 2r real fields. 

The DT methods would be applied following similar steps as in section 3 and subsection 3.1, for constant 
and free field NVBC's, respectively. In particular, in the constant NVBC the form of the relationships 
(3.25) and (3.26)-(3.28) will remain the same, except that i = 1, 2, r. In the case of free field NVBC the 
relationships (5.47)-(5.48) would be satisfied with i = 1,2, ...,r. 

The VBC and the bright solitons will be associated to the vertex operators Fj,Gj (see (C.l)) as in 
section 4. In this case, one requires pf = 0, (i = 1,2, ...,r) in (3.25). The dark solitons, as in section 
5, will require the vertex operator of type W q (k,pf) (j = 1,2, ...,r), the analog of the operator in (C.7) 
incorporating additional terms. Finally, the mixed boundary conditions and the dark-bright solitons will 
emerge by extending the discussion in section (6). In the case of the vector 1-soliton solution it is possible to 
form the combination (m, r — m), m ^number of dark components, r — m =number of bright components. 
So, the vertex operators analog to r± ( g 2 (/c ± , pf) in (C. 3) will be associated to the roots ±/3j, (±f3j TPi) (i = 
1,2,...,to) such that T±p j {k ± ,pf) (j = m + 1, ra + 2, r) . 

8 Discussion 

We have considered soliton type solutions of the AKNS model supported by the various boundary condi- 
tions (3.2)-(3.5): vanishing, (constant) non- vanishing and mixed vanishing-nonvanishing boundary conditions 
related to bright, dark and bright-dark soliton solutions, respectively, by applying the DT approach as pre- 
sented in [15]. The set of solutions of the AKNS r system (2.6)-(2.7) is much larger than the solutions of the 
r-CNLS system (2.11). A subset of solutions of the AKNS r system, (2.6)-(2.7) for r = 2 and (5.4)-(5.5) for 
r = 1 , respectively, solve the scalar NLS (5.7) and 2— CNLS system (2.11), under relevant complexifications. 

Moreover, the free field boundary condition (3.29) for dark solitons is considered in the context of a 
modified DT approach associated to the dressing group [29], and the general N-dark-dark soliton solutions 
of the AKNS2 system have been derived. These soliton components are not proportional to each other and 
thus they do not reduce to the AKNSi solitons, in this sense they are not degenerate. We showed that these 
solitons under convenient complexifications reduce to the general N-dark-dark solitons derived previously in 
the literature for the CNLS model [4, 14]. In addition, we have shown that two— dark— dark— soliton bound 
states exist in the si (3)— AKNS system, and three— and higher— dark— dark— soliton bound states can not 
exist. These results hold for any value of the index q. In the case of reduced 2— CNLS when focusing and 
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defocusing nonlinearities are mixed, this result corresponds to 2-dark-dark soliton stationary bound state 
[4]- 

In the mixed constant boundary conditions we derived the dark-bright solitons of the sZ(3) AKNS model. 
These solitons under the complcxification (6.1) reduce to the solitons of the 2-CNLS model (6.2) which will be 
useful in order to investigate dark-bright solitons appearing in an inhomogeneous two-species Bose-Einstein 
condensate [26]. 

The relevant steps toward the AKNS r (r > 3) extension were briefly discussed in the framework of the 
DT methods. In particular, the vertex operator calculations can be extended in a direct way following the 
same steps as in the appendices B and C and the sl(n) highest weight representation [42]. 

Another point we should highlight relies upon the possible relevance of the CNLS tau functions to its 
higher-order generalizations. We expect that the tau functions of the higher-order CNLS generalization are 
related somehow to the basic tau functions of the usual CNLS equations. This fact is observed for example 
in the case of the coupled scalar NLS+ derivative-NLS system in which the coupled system possesses a 
composed tau function depending on the basic scalar NLS tau functions [41]. 



9 Acknowledgments 

The authors thank the referees for relevant comments and suggestions and M. Zambrano for discussions. HB 
has been partially supported by CNPq. AO A thanks the support of the brazilian CAPES. 



A si (2) matrix elements 

The commutation relations for the sl(2) affine KacMoody algebra elements are 

[ H (m) fH (.n)] = 2mS m+rh0 C, 



H (m) ^ E (n) 
■ E (rn) ^ E (n) 



d , ri m ) 



= ±2E. 



(m+n) 



= ("*+") + md m+n . C; 

= mTi m ); = {H^ , E^} 



(A.l) 
(A.2) 
(A.3) 
(A.4) 



The central extension ensures highest weight representations (h.w.r.) of the affine algebra (sec e.g. [15]). 
So, in the h.w.r. {|A >, |Ai >} one has the following relationships 



4 0) K> = o 



4 m) K> 



0, TO > 

#(™)|A a > = 0, m>0; 
ff (0) |A a > = MA a >, 
C|A a > = |A a > 



(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 
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where a = 0, 1. The adjoint relations (^"V = E ( T m) , {H {m ^ = H^ m '> allow one to know their actions 
on the < A a |. Next, consider the vertex operators 

v«h,p) = E i(7 2 - p 2 r n/2 Ki" - - _ eq ^ 2 ~ 2)1/2 4 W) + 



n— — oo 



-E ( 

1 + e q {^-py/* - 



+ ^( Z ^ 2 JL ) 1/2 KoC}; 9 = 1,2; (A.10) 
where e q = (— and p 2 = f> + f>~ ■ The vertex operator V q (j, p) satisfies 

[e 1 ,V*]=* r V*, [s 2 ,V q ] = (-l) q - 1 2 1 { 1 2 -p 2 ) 1 l 2 V q , q=l,2, (A.ll) 

where 

fx = ff« + /5+4 0) + p-i?i 0) , e 2 = + p+eW + p-E W . (A.12) 

The following matrix elements can be computed using the properties (A.5)-(A.9) 

h 2 - p 2 ) 1/2 

(AolV'IAo), = e q K -l p ->—, (A.13) 

7 

<A |£«^|A ) = T - X ~ 2W2 (J 2 -P 2 ) 1/2 - (A.14) 
7T &q\T - P 2 Y' 2 

The matrix element (A | V q (ji, p)V q (j 2 , p) |A D ) can be computed by developing the products and keeping 
only non-trivial terms, then one makes use of the commutation rules to change the order, and eventually 
to get some central terms C . The double sum can be simplified to a single sum and each term can be 

substituted by power series like J2™=o x " = (T^-y , J2^=i x " = Ji^x) ' an ^ S^Li nx " = (i-x) 2 • one can 
get 

(A |^( 7l ,p)^( 72 ,p)|A ) = [2 + 2 g(^^)] ^^L- + 

[ X(71 ' 72) ^ +1]; , = 1,2 (A.15) 

7 1 7 2 

where tf ( 7l>7 >) = - ^IZ^IZZ^ (A . 16) 

Vif - p 2 p 

In order to prove the nilpotency property of the vertex operator V q ( 7 i , p) , when evaluated within the 
state |Ao >, it is convenient to write (A.15) in the following Laurent series expansion 

^2 



1_ 10 72 (71-72) 15(6 72 2 + ;5 2 )( 7 i- 72 ) 2 
4! (7 2 2 -/5 2 ) 5! + hl-p 2 ) 2 6! 
420 72 (2 7 2 + /) 2 ) ( 7 i- 72 ) 3 



( 7 2 - p 2 f 7! 

From this, it is clear that 



+ ... 



(A.17) 



lim (A |^(7i,/5)^(72,/5)|A ) 0, q = 1,2. (A.18) 

7l->72 
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B The affine Kac-Moody algebra sl 3 (C) 

In the following we provide some results about the affine Kac-Moody algebra Q = sh(C) relevant to our 
discussions above. We follow closely [42, 18]. The elements of the sl 3 (C) Lie algebra are all 3 x 3 complex 
matrices with zero trace. Consider the corresponding root system A = {±ai, ±a 2 , ±"3}, such that the three 
positive roots are a*i, i = 1, 2, 3, with a a , a = 1, 2, being the simple roots and «3 = ot\ + a 2 . We choose a 
standard basis for the Cartan subalgebra % such that 

/ 1 n n \ / ^ 

1 
0-1 



#1 = 



1 
0-10 




H 2 



\ 



(B.l) 



and the generators of the root subspaces corresponding to the positive roots are chosen as 



( n 1 n \ / \ / 1 \ 



E 



+Ctl 



1 





E. 



+ C2 




1 




E. 



+ U3 



1 





(B.2) 



For negative roots one has E- a = (E +a ) T . The invariant bilinear form on sh{C), (x\y) — tv(xy);x,y € 

sls(C) induces a nondegenerate bilinear form on H* which we also denote by (• | •). This definition allows 
one to write 

(ai|ai)=2, (a 2 |a 2 ) = 2, (ai|a 2 ) = -1. (B.3) 

On the other hand, the generators = {H[ m \ H^, E^}, where m <E 7L and a e A, together with 
the central C and the 'derivation' operator D {[D, T^] — mT^>) form a basis for sZ 3 (C). These generators 
satisfy the commutation relations 



D-(m) rr{n) 
n a ! n b 

H (m) E (n) ' 
E (m) £.(«) 

"a ' —"a 

R (m) R (n) 

«3 ' — OL3 

p{m) p(n) 
"ai ' 0:2 

p;(m) r(n) 
«3 ' -ai 

E (m) tp(«) 

[£>, C 



m (a a |ab) C <5 
±(a a |a l )^ 



J ±OLi 



_ K (m+n) 
R("»+n) 

0, 



(m+n) 



+ m C 5 m+nj o? 



(B.4) 
(B.5) 
(B.6) 
(B.7) 
(B.8) 
(B.9) 
(B.10) 
(B.11) 



where a, 6 = 1,2, i = 1, 2, 3 and m, n e 2. The remaining non-vanishing commutation relations are obtained 
by using the relation [ E^ , £?i n) 1 = - \ si~ m) , 
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In this paper we use the homogeneous Z-gradation of sl 3 (C) which is defined by the grading operator 
D, such that 

sh{C) = Gm, [Gm ,Qn]C Gm+n, (B.12) 

me ~2- 

Where Q m = {x e sl 3 (C) \ [D,x]=mx; m <G Z}. 
The subspace Go is a subalgebra of sl 3 (C) given by 

g = CH 1 @CH 2 @CE^ ®CC®CD (B.13) 

and for the subspaces Gm {m ^ 0) we have 



Gm = c H[ m) e c H { 2 m) e c £<™> e c e c e^ ® c e { Z{ e c £^ } 2 e c E_ a 



(m) 

(B.14) 



We use in the paper the fundamental highest weight representation | A ) , satisfying 

ff(°)|Ao)=0, £?(°)|Ao) = 0, C|Ao) = |A ) (B.15) 

for a, b = 1, 2, and a e A. Such state is annihilated by all positive grade subspaces 

g m |A )=0, rn>0, (B.16) 

and all the representation space is spanned by the states obtained by acting on |Ao > with negative grade 
generators. This representation space can be supplied with a scalar product such that one has 

(H™)i = H~ m , (2^)t = EZl\ (B.17) 
C+ =C, Dt = £>. (B.18) 

It follows from (B.13) and (B.14) that 

{Qm) ] = G-m, (B.19) 

and therefore 

<A o |£L m =0, m>0. (B.20) 
In addition to the subalgebra g it is also convenient to consider two additional subalgebras 

G<0 = G m , G>0 = G-m- (B.21) 

m>0 m>0 

These subalgebras and the corresponding Lie groups play important role in the DT method. 
The next relationships are useful in the AKNS r (r = 2) model construction. The special element E^ in 
the basis presented above can be written as 

= l{H[ l) +2H^), [D,E®] = IE®. (B.22) 
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Then the matrix becomes 



£<°) = - 
3 



( I 



V 



I— 1 











1 











-2 



/ 



The roots entering in the AKNS2 construction are 

Pi = 013 = ai + a 2 ; j3 2 = a 2 . 
Moreover, the following commutation relations hold 



[EW,HW]=l52aC6 l+m , , a =1,2, 



1 ,-,(Z+m) • _ 1 9 

= 0, 

= -|-_£;( m +™) 



if. 



(to) ,-,(71) 



2 i ^±/3 



l_ p( m + ™) 



^1 ' tj ±h 



(to) 5,(71) 
2 1 



= ±2£. 



(m+n) 
±02 



#1 



(to) 5,(71) 



^2 1 Ej ±ft^fa 



±2S 



(m+n) 



E (m+n) 



(B.23) 



(B.24) 

(B.25) 
(B.26) 
(B.27) 
(B.28) 
(B.29) 
(B.30) 
(B.31) 



C s/(3) matrix elements 

Consider the vertex operators associated to bright soliton solutions 



+00 



+00 



(C.l) 



It can be shown that they are nilpotent, i.e. F? = 0,G 2 = 0. The matrix element (X \FjGk \X a ) can 
be computed by developing the products and keeping only non-trivial terms, then one makes use of the 
commutation rules to get the central term C. The double sum can be simplified to a single sum, which 



provide the power series Y^=i nx ™ = (i-s) s • one nas 



(A | FjGk |A D ) 



Pfe 



- Pk) 

Let us consider the deformation of the vertex operators F 2 , G 2 as 



+00 + 
,10 



(C.2) 



(C.3) 



It is a direct computation to show the nilpotency of these operators, i.e. ^±^ 2 = 0. Similar computations 
to the one in (C.2) provide the following matrix element 

w + w~k + k~ 



(X o \T 02 (k + )T^ 2 (k-)\X o ) = 



(k+k~ +ptpi)(k+ -k~) 2 



(C.4) 
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Consider the vertex operator analog to the one in (A. 10) 

1_ 

X-e q (A2-p2)i/2^ 



71 — — OO 



where e q = (— 1) 9_1 and p\ = 4p^p^. The next matrix element computation follows similar steps to the one 
performed to arrive at (A. 15), except that one must take into account the sl(3) commutation rules. So, one 
has 

(A |^ i (A 1 ,po)^ 1 (A 2 ,p )|A ) = \{[2 + 2K{\ U A 2 )] ^ ^^l 2 )? + 

[ K{X li 2)P K l]}; 4=1,2, (C.6) 

where Ko&ndK are given in (A. 16). Since this two-point function, except for an overall constant factor, is 
similar to the one in (A. 15) one can use the relationships (A.17)-(A.18) to show that the operator (Ai, po) 
is nilpotent. 

The vertex operator generating the dark-dark soliton solution becomes 

00 2 2 

w"(*,p£ a ) = £ {(fc 2 -4E^r)-" /2 [^"[^^ ri) + ^ ri) + E< 9 )< ) + 

n— — 00 2—1 i—1 

f e - E {n) +e+E {n) +e-E (n) 1 1 c ( fc2 ~ 4 EL P+/Q 1/2 5 g i 

i (?) -ft + 12^/31-/32 + e 12 £/ /3 2 -/3 1 J + e 9 2 jfc °n,0W 

(C.7) 

+ _ tp? n _ 1. o _ 1 Pi>r . „± _ 1 P1P2 . 



(-ir 1 



Notice that the vertex operator (C.7) reduces to the one in (C.5) in the limit p^ — > 0. The nilpotent 
property of this vertex operator can be verified as follows 

(Ao| W\k u p± 2 )W\k 2 ,p± 2 ) |A ) = (^l)^l ^ + (»2 - *i) + •••) (C8) 

where xi = \Jk\— 45, x 2 = — 4S*, S* = ^ • Pj^pJ ■ In the limit a; 2 — > £1 (or fc 2 — > fci) the r. h. s. of eq. 
(C.8) vanishes. 
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